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DISCLAIMER 



Although each program has been tested by its contributor, no 
warranty, express or implied, is made by the contributor or 
1620 USERS Group, as to the accuracy and functioning of the 
program and related program material, nor shall the factor 
distribution constitute any such warranty, and no responsibility 
is assumed by the contributor or 1620 USERS Group, in con- 
nection therewith. 



Modifications or .revisions to this program, as they occur, 
will be announced in the appropriate Catalog of Programs 
for IBM Data Processing Systems. When such an announce- 
ment occurs, users should order a complete new program 
from the Program Information Department. 
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Are the Sense Switch options adequately described (if applicable)? Yes No_ 

Are the mnemonic labels identified or sufficiently understandable? Yes No^ 
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ABSTRACT 
LINEAR PROGRAMMING II 
by 

F . tf. flood 



DESCRIPTION 

This Linear Programming code is essentially the same as the previous 
four-phase code, the difference being the addition of an output phase and 
the reduction in the number of typeouts in the other phases. This writeup 
assumes familiarity with the four-phase writeup, the difference being what 
is included here. 



DISCUSSION 

The revision of the four-phase program to the five-phase program was 
initiated by the need for extra matrix capacity. Error typeouts were 
eliminated and only the pause now exists. The output portion of the program, 
after the optimum solution is reached, was removed from the Solution Phase 
and from the Right Hand Side Changer Phase and made an Individual Phase. 
Therefore, to obtain a solution it is necessary to lodd the Output Phase after 
the final stop has been encountered in the Right Hand Side Changer or the 
Solution Phase. 

The MAR display number as described in the following sections are all 
based on FORTRAN with auto-divide and a lowered origin. If the programs are 
compiled with FORTRAN decks that do not comply with the above, these numbers 
will differ from those reported below and the correct number must be established 
from the print off when the program is compiled. 

DATA LOADER PHASE 

Input 

Same as L.P.I. 
Output 

Same as L.P.I. 
Sense Switch Settings 

On Off 
SS 1 For restart problems For initial start problem^ 

SS 2,3,4 Not used 
Error Stops 
None 
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SOLUTION PHASE 
Input 

None 
Output 

Same as L.P. 1 except no headings and no final output. 
Sense Switch Settings 

On Off 



SS 4 



SS 1,2,3 
Error Stops 
MAR Display 
9039 
12011 
12023 
COST CHANGER PHASE 
Input 

Same as L.P, 1 
Output 

Same as L.P. 1 
Sense Switch Settings 

None used 
Error Stops 

None 

RIGHT HAND SIDE CHANGER PHASE 
Input 

Same as L.Po 1 



Punch output for restart at a 
later time 

Not used 



Continue with Simplex 
Method 



Meaning 

Problem as stated is unbounded 
Problem as stated is infeasible 
Optimum solution has been reached 
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RIGHT HAND SIDE CHANGER PHASE (Cont'd) 
Output 

Same as L.P. 1 with the exception of error typeouts and the 
final solution. 

Sense Switch Settings 

None used 
Error Stops 

MAR Display Meaning 

10651 Problem is infeasible 

12575 Optimum feasible solution has been reached 

An alternate method of input is available without any program change. 
The original method of input was to punch both the old right hand side element 
and the new right hand side element into the change card. 

It is possible to record only the difference between the elements in card 
columns 9-18 and leave the rest of the card blank. 

The card should then be punched as follows: 

cc 

1-4 l) 0000 if the unit vector associated with tne change is artificial 

2) 0001 if the unit vector associated with the change is slack 

5-8 Column number of unit vector associated with the change 

9-18 Old right hand minus the new right hand side with a decimal 

point somewhere in the 10 digit field. Care must be taken 
with the sign of the resultant of the subtraction. To 
increase the right hand side element, the number must be 
negative and to decrease it must be positive. 

OUTPUT PHASE 
Input 

None 
Output 

Solutions shadow costs, final matrix punch, and solution punch 
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Sense Switch Settings 

On 

SS 1 Print final solution 

SS 2 Punch final matrix 

SS 3 Print Shadow Costs 

SS 4 Punch final solution 

In all the above cases the final 
Error Stops 
None 



Off 

Do not print final solution 
Do not punch final matrix 
Do not print Shadow Costs 
Do not punch final solution 

solution is printed on the typewriter. 



Mathematician 
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g '-:.co 
c . eo 

O.'.'OOO 
O - O2o 

0150 
":i 26 
: )193 
00246 
03294 
03306 

g'^10 

03^66 
03350 
03353 
00632 
Oo730 
GO/36 

c307o 

08913 
03926' 
03 ; Uf6 
0:5970 
09006 
C9074 
09154 
09166 
09282 

03493 

C9534 
095^2 
09566 
09578 



DIMENSION A(44,99),'./(44),L(44),I I ( 1 ) JJ ( 1 ) , I I I ( 1 ) ,X( 1 ) 

1 F0-MATO4 f^ 1) ' KJ(1) ' KKK(1) ' IJ(1) ' K ^ 1) » , ( 1 >' J ^)» XM,N (1) 

2 FORMAT ( l4,*i4,F20.8) 

4 FORMAT ( 1 5HL0AD NEXT PHASE ) 
108 PEAD 1 , t I , J J 
I I 1 = 1 1+1 
DO 10 ( = 1,11 I 

V/(!)=0.0 - - 

L(l)=0 

DO 10 J = 1 ,JJ - - ~ 

10 A( I ,J)=0.0 

11 READ 2,! ,J,X - 
1F( I ) 1 6 , 1 6 , 1 8 

18 A( I ,J)=X - -. - 
GO TO 11 

16 J=JJ • - - 

17 READ 2,i ,JK,X 

iF{ 1)20,20,19 - -~ — 

19 A( I ,JJ) = X 

L( t )=JK --• — - - - 

GO TO 17 

20 I F ( SENSE SWITCH 1)40,22 

22 1=1 

23 1=1+1 " - " 

IF( l-i I I ) 24 ,40 ,40 

24 IF(L(J ))23,25 ,23 ~ 

25 DO 27 J = 1 ,JJ 

IF(A(I ,J))26,27,26 - 

26 A( I 1 I ,J)=A(I ! I ,J)-A(i ,J) 

27 CONTI NUE - 

GO TO 23 

40 PR I NT 4 - - 

PAUSE 

END - " - 
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SOLUTION PHASE 



0-000 DIMENSION A ( 44 , S'S ) , ■ ■' ( 44 ) , L ( 44 ) , I I ( 1 ) , J J ( 1 ) , I I I ( 1 ) ,X( 1 ) 

03000 DIMEHSIOM J K ( 1 ) , K J ( 1 ) , KKK ( 1 ) , IJ(1) ,K(1 ), I (1) , J ( 1 ) , X M I N ( 1 ) 

O.iOOO 105 FORMAT (I4,F15.2,I4) 

03032 !<KK=0 

03056 K— I I I - 

03030 44 J=0 

03104 W(K)=0.0 

03152 L(K)=0 

03200 42 J=J+1 " " 

03236 IF{J-JJ)41 ,45,45 

03304 41 IF(A(K,J)143,42,42 

03420 43 iF(W(K)-A(K,J))42,42,47 

03572 47 V/(K)=A(K,J) 

00630 L(K)=J 

03723 GO TO 42 

03736 45 IF(L(K))46,62,46 

03316 46 KJ=L(K) 

03364 DO 120 I =2, 1 I 

03376 I F( A( I ,KJ ) ) 1 20 , 1 20 , 1 2 1 

03992 120 CONTINUE 

09023 PAUSE 

09040 121 1=1 

09064 ~ JK=0 ' - 

09038 50 1=1+1 

09124 IF( l-l I )52,52,56 

09192 52 IF(A(i ,KJ»50,50,51 

09308 51 X=A( I ,JJ)/A(I ,KJ) - 

09464 IF(JK)55,53,55 

09520 "-55 IF(X-XMtN)53, 50,50 - 

09538 53 XMUi=X 

09612 JK=I ■ -■- 

09636 GO TO 50 

09644 56 X=A(JK,KJ) - - 

09728 L(JK)=KJ 

09776 DO 57 ! =1 , ! I I 

09788 57 V/(l )=A(I ,KJ) 

09932 IJ=JK-1 

09968 DO 59 1 = 1 ,U 

09980 DO 59 J = 1 ,JJ 

09992 1F(A(JK,J))58,59,58 

10108 58 ! F(V/( I ) ) 580, 59, 580 

10188 580 A(l ,J)=A(I ,J)-W(I)*(A(JK,J)/X) 

10464 59 CONTINUE 

10536 IJ=JK+1 

10572 DO 61 l = !J,l I I 

:3534 DO 61 J=1,JJ 

10596 IF(A(JK,J))60,61 ,60 

T0712 60 IF(W(I ))600,61 ,600 

10792 600 AO ,J)=A( I ,J)-WO )*(A(JK,J)/X) 

11068 61 CONTINUE 



7 

SOLUTION PHAS-: (COMT) 



JM 140 DO 205 J = 1 ,JJ 

11152 205 A(JK,J)=A(JK,J)/X 

H344 KKK=KKK+1 

11380 PRINT 105,KKK,A(K,JJ) ,L(JK) 

11512 IF(SENSE SWITCH 4)70,44 

11532 62 IF(K-1)70,70,63 

11 600 63 I J=JJ-1 

11636 DO 65 J=1 ,U 

Tl648 !F(A(K,J)-.0001)65,65,66 

11776 65 CONTINUE 

11312 DO 130 J=1,JJ 

Tl824 130 A(lll,J)=0.0 

11944 K=1 

T 1 963 KKK=0 

T1992 GO TO 44 

12000 66 PAUSE 

12012 70 PAUSE 

12024 END 
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\ GHT HAND S i DE CHANGER 



o;:ooo 

Ow-000 

coooo 
00033 

0312Q 
00156 
03244 
03304 
03328 
00352 
00376 
03460 
08516 
00572 

; ors596 
03604 
011628 
08688 
03748 
O876O 
09024 
09032 
09056 " 
09080 
09128 
09152 
09188 ~ 
09256 
09372" 
09500 
09584 
09632 
09640 
09720 
09768 
09792 
09816 
09364 
09900 
09968 
T0084 
T0 188 
10332 
70412 
T0480 
T0504 

"10552: 
10560 
10640 
10652 



0!. HEMS I ON A( 44, 99) ,V/(44) ,1.(44) , I I (1 ) ,JJ( 1 ) , I I I (1 ) ,X{ 1 ) 
D I MENS ! ON J K{ 1 ) , KJ { 1 ) ,K'\K( 1 ) , !J ( 1 ) ,K( 1 ) , I ( 1 ) , J ( 1 ) , XM I N ( 1 ) 
1 FORMAT ( 14, l4,F10.5,n0.5) 

3 FORMAT (33 H ITER FUNCTIONAL. VAR. IN) 

100 F0RfiAT(F20.5,F20.5) 



101 FORMAT ( 32 H OLD PHS 

104 FORMAT ( i 4 , F 1 5 . 2 , 1 0X , I 4 ) 

3 PRINT 101 
K=1 
KKK=0 

4 READ 1 .KJ.JK.Wd) ,W( ill) 
IF(KJ)81 ,80,81 

00 IF(JK)82,11 ,82 

81 IJ = 1 

GO TO 5 

82 IJ=2 

5 PRINT 100,W(1),W(I II) 
W(1)=W(M I )-W(1) 

DO 10 l=IJ,l I 

10 A(l ,JJ)=A(I ,JJ)+W(1)*A(f ,JK) 
GO TO 4 

11 PRINT 3 

"15 XMIN=0.0 " 
L(K)=0 
1=1 

18 [=1+1 

IFff-l 1)16,16,20 
16 IF(A(t,JJ))17,fl3,l8 
"17 "IF(A(l,JJ)-XMIft)t9, 18,18 — 

19 XMIN-A(I ,JJ) 
L(K)=f 

GO TO 18 

20 IF(L(K))21 ,70,21 

21 JK=L(lO 

'"'*' j=o - " " ' 

XMIN=0.0 
L(K)-0 
26 J=J+1 

1F(J-JJ)23,24,24 

23 IF(A(JK,J))25,26,26 
25 IF(A(1,J))26,26,7 

7 X=A(1 ,J)/A(JK,J) 
1F(L(K))6,27,6 

6 IF(X-XHIM)26,26,27 
" 27 XMIN-X 

L(K)»J 
GO TO 26 

24 IF(L(K))31,30,31 

30 PAUSE 

31 KJ=L(K) 



NEW RHS) 



JO 



6 



o 



RIGHT HAND SiDE CHANGER (CONT) 



10700 
TC748 
TC832 
70344 
To 938 
11024 
11036 
11048 
31164 
1 1 244 
11520 
11592 
11628 
11640 
T1652 
71763 
11843 
12124 
12196 
12208 
12400 
T2436 
T2556 
12564 
32576 
12584 



L(JK)=KJ 
X=A(JK,KJ) 
DO 33 1=1,111 
1/(1 )=A(I ,KJ) 
!J»JK-1 
DO 35 1=1 ,U 
DO 35 J==1 ,JJ 
IF(A(JK,J))34,35,34 



34 IF(W(I))37,35,37 

37 A(l ,J)-A(I ,J)-W(I)*(A(JK,J)/X) 

35 CONTINUE 
U=JK+1 

'-DO 40 l.«IJ,l II 
DO 40 J-1.JJ 
IF(A(JK,J))38,40,38 

38 IF(W(I))39,40,39 

39 A(l ,J)=A(I ,J)^W(I)*(A(JK,J)/Xr 

40 CONTINUE 

DO 50 J=1 ,JJ 

50 A(JK,J)=A(JK,J)/X 
KKK=KKK+1 

PRINT 104,KKK,A(1,JJ) ,L(JK) 
GO TO 15 
70 PAUSE 

GO TO 8 " - " 
END 
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COST CHAFER 



00000 D I MEMS I ON A(44 ,99) ,W(44) ,L(44) , I I ( 1 ) ,J J( 1 ) , I I I (1),X(1) 

08000 DIMENSION JK(1) ,KJ(1 ■ ,KKK( 1 ), U ( 1 ) ,K( 1) , I ( 1 ), J ( 1 ), XMI N( 1 ) 

03000 1 FORMAT ( l4,F10.5»F10o) 

03032 2 FORMAT ( 13HSTI LL OPTIMAL) 

0G082 3 FORMAT ( 10HFUNCTI ONAL ,F20.8) 

00132 4 FORMAT (18HL0AD SOLUTION DECK) 

08192 100 FORMAT ( 25 HVAR. OLD COST NEW COST ) 

08266 9 PRINT 100 

08290 10 READ 1 ,JK,V/(1),W( III) ~ 

08362 L(1)=JK 

08386 IF(JK)11 ,50,11 

08442 11 PRINT 1,JK.W(1),W(II I) 

08514 W( 1 )=W( I I I ) -W( 1 ) " 

08574 A(1 ,JK)=A(1 ,JK)+W(1) 

08706 DO 15 1=2, 1 I - : " 

08718 IF(L(I)-L(1))15,16,15 

08810 ' 15 CONTINUE ; ■- " " 

08846 GO TO 10 

68854 16 DO 20 J=1,JJ ~ ~ 

68866 IF(A(I ,J))19,20,19 

08982 19 A(1,J)=A(1,J)-W(1)*A(I,J) 

09198 20 CONTINUE 

' 69234 GO TO 10 

09242 50 I J=JJ-1 

09278 '" , - DO 55 J=l , I J 

09290 IF(A(1,J))60,55,55 

09394 " 55 CONTINUE ~ 

09430 PRINT 2 

09454 PRINT 3,A(1,JJ)~" " " 

09526 PAUSE 

09538 GO TO 9 ~"~ 

09546 60 PRINT 4 

09570 PAUSE 

09582 GO TO 9 

09590 " END 



If 

OUTPUT ; J .-vS'ii 



c.';ooo dimension A{44,99) ,v/<;44) ,l(44) , i i I (1),x(i) 

C:<CCO D ! MENS I ON J!<C 1 ) ,KJ ( 1 } ,KK'<( 1 ) , U ( 1 ) ,K( 1) , I ( 1 ) , J ( 1 ) , XMI N( 1 ) 

COOCO 2 FORMAT( 14, I4,F20.3) 

C1C32 5 FORMAT ( I if ,F20.3) 

£2-060 7 FORMAT ( 23HVARI ABLE VALUE) 

0":150 8 F0RMAT(10HFUNCTI0NAL,F20.3) 

03l£0 9 FORMAT ( 23 HVAR I ABLE SHAD . COST) 

03 250 100 FORMAT ( 1 H ) 

03?7° I F( SENSE SWITCH 1)70,92 

CT;296 70 PRINT 3,A(1,JJ) 

PRINT 7 

03392 DO 71 1=2, 1 I 

0-3404 71 PRINT 5,L(I),A(I ,JJ) 

03560 92 !F(SENSE SWITCH 3)72,75 

03530 72 PRINT 9 

03604 !J=JJ-1 

03640 DO Ik J = 1 ,U 

08652 IF(A(1 ,J))73, 74.73 

OG756 73 PRINT 5,J,A(1,J) 

0R340 74 CONTINUE 

03376 75 i F( SENSE SWITCH 2)76,80 

03396 76 IJ=JJ-1 

03932 DO 77 J = 1 ,U 

03944 DO 77 1=1 J M 

60956 IF(A(I ,J) 573,77,73 

09072 78 PUNCH 2,1 ,J ,A( I ,J) 

09130 77 CONTINUE 

09252 PUNCH 100 

09276 91 DO 90 I =1 , 1 M 

09238 90 PUNCH 2 , I ,L( I ) ,A( I , JJ ) 

09456 31 PAUSE 

09463 80 IFCSENSE SWITCH 4)91,31 

09438 END 
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